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sibility to discriminate models of the charm hadroproduction including the low-x QCD be- 
haviour of hadronic cross-sections relevant at very high energies. 
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INTRODUCTION 
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The cosmic ray (CR) spectrum is characterized by a sharply falling power law behaviour, 
£;-(7+i)j2|. The spectrum gets more steeper around 10^ GeV with the spectral index 7 changing 
from 1.7 to 2.1 - this region is called the knee. Around ~ 5 x 10^ GeV, one observes a flattening 
of the spectrum, with the spectral index 7 falling between 1.4 and 1.7. This is the so called ankle. 
These two breaks of the primary spectrum are still open questions of CR physics. The region 
beyond the ankle is the regime of ultra high energy cosmic rays. There is not much data available 
in that region and no clear consensus exists on the composition or the particle content in this region 
[^. It is generally believed that the change in the slope around the knee is astrophysical in nature 
rather than any specific change in hadronic properties and/or interactions j, y]. An overview of 
hadronic interactions and cosmic rays can be found in [4']. 

The atmospheric muon flux originated from decays of pions and kaons is commonly called the 
conventional muon flux. There is expected rather sharp reduction of the conventional muon flux 
above a few TeV [3] due to the increasing decay lengths and decreasing interaction lengths of 
pions and kaons. Therefore, at very high energies the bulk of muons is expected to arise from 
the semileptonic decay modes of heavy shortlived hadrons, predominantly the charmed ones. This 
component is called the prompt muons. It is known that the prompt muon flux is only about 
10% smaller than the prompt flux at the surface of the earth. Therefore, measurement of 
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the atmospheric prompt muon (APM) flux at high energies wih ensure a normaUzation for the 
atmospheric prompt neutrino (APN) flux, and a direct comparison of the two is both desirable and 
necessary. This study is necessary because the atmospheric neutrino flux is unavoidable background 
to VHE neutrino experiments. 

There are sizeable uncertainties in theoretical predictions for the prompt lepton fluxes (see P, Q] 
for review). The reason is mainly due to the vastly different choices for the charm production 
cross-section - perturbative QCD (pQCD) with a K factor [8], next-to-leading order (NLO) pQCD 
phenomenological nonperturbative approach, such as the recombination quark-parton model 



(RQPM) or quark-gluon string model (QGSM) [6|]. The experimental situation is not very precise 
either at this stage. Various experiments provide upper limits on the APM fluxes in the energy 
range of interest, which allow a large variation in the prompt fluxes. One can therefore expect that 
better measurements of high-energy muon fluxes can play a definitive role in selecting the charm 
production models, and thereby, also providing invaluable information about parton densities at 
such low-x and high energy values. Another related source of large theoretical uncertainties is 
strong dependence of the hadronic cross-sections on the renormalization and factorization scales. 
This is partly related to the naive extrapolation of parton distribution functions to very different 
energy and x-values. For the case of conventional fluxes originating from the pions and kaons, 
these issues are in much better control and therefore the predictions stand on a sound footing. 

Earlier authors of Ref. [l^, [l^ explored the possibility of utilizing the high energy prompt muon 
flux(es) in order to investigate whether the general expectations expressed above can in practice 
help in selecting the charm production model/parameterization and also the importance of the 
heavy composition of cosmic rays above knee. They chose some of the models often used and 
compare the predictions, incorporating the saturation model of Golec-Biernat and Wuthsoff 
However while esimating their event rates of muons in a 50 kT Iron detector like INO one [l^ they 
did not consider the angular dependence of the muon fluxes at rock depth. Angular dependence 
of muon flux due to surrounding rock is really important for correct estimation of the muon event 
rate inside such a detector. In this work we calculate the high-enery AM flux, conventional as well 
as prompt, at INO rock depth taking into account the distortion in the surface muon zenith- angle 
distribution due to specific topography of the INO site. 

It is therefore quite clear from all these models that the lepton fluxes at the end are strongly 
sensitive to the charm production cross section. Till the knee, the cosmic ray flux and composition 
is rather established and therefore, the only source of large error is the charm cross section. This 
therefore gives us a unique possibility to gain information about heavy quark production mechanism 
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at high energies and low x. 

II. SURFACE ATMOSPHERIC MUON FLUX AND THE CALCULATION TECHNIQUE 

A. Topography of PUSHEP site 

The slant depth X depends on the topography of the rock surrounding the INO detector. 
PUSHEP is the selected site for this purpose. One can assume a constant depth which is equal to 
the vertical depth just above the cavern. The vertical depth of PUSHEP site is 1.3 km of rock. 
Another assumption is that of a triangle topogarphy. In this case the slant depth for given zenith 
angle 6 is calculated as 

ho 



X{9) 



(1) 



cos 6* + (ho/lo) sin( 

where /iq = 1-3 km is the vertical depth, /q = 2.1 km is the half-length of the approach tunnel and 
tanuj = Hq/Iq is the slope of the mountain. The triangle nature of site and the slant depth X{9) 
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FIG. 1: Geometry and slant depth of PUSHEP site. 

are shown in Fig. [TJ For the rock density g we adopt here value 2.72 g/cm^. The column depth 
h{9) related to the slant depth, h{6) = qX{9), varies between /imin — 3.0 • 10^ g-cm~^ (3 km w. e.) 
that corresponds to cos^m — 0.85 and ^max — 5.71 • 10^ g-cm~^ (5.71 km w. e.) near horizontal. 
Near vertical direction column depth is about 3.54 km w. e. 



B. Parameterization of the conventional muon spectrum at sea level 



The surface muon flux is rather well measured up to TeV and can be described by different 
analytical formulae taking into account the zenith-angle dependence. Here we list some of them 



which were used in present calculations. First of all we use Gaisser's muon flux parameterization 
[l6 | (in inits of cm^^s^^sr^^GeV^^) 



-2.7 



+ 



0.054 



1 + l.l(£;^/115GeV)cos6' 1 + l.l(£;^./850GeV) cos 6^ 
For our purpose we work with a modified muon flux formula obtained by Tang et al. 



(2) 
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Next parameterization of the conventinal muon flux we use here is that by Bugaev et al. 
for vertical direction: 



1^ 



^I'^iPt^,^! = , cm-2s~isr'i(GeV/c] 



-1 



(3) 



where z = logio(p^t/l GeV/c). Values of parameters in Eq. ([3|) are listed in Table H] for different 
momentum ranges. The muon energy spectrum is 4>li^^{E,9) = {E^/pfj) (j)li^^{p^,6). For inclined 

TABLE I: Parameters in Eq. ^ for the vertical energy spectrum of conventional muons at sea level. 



Momentum range, GeV/c C, (cm^s sr GeV/c) ^ 70 



71 



72 



73 



1.0 - 927.65 
927.65 - 1587.8 
1587.8 - 4.1625 x 10^ 
> 4.1625 X 10^ 



2.950 X 10^3 
1.781 X 10-2 
14.35 
10^ 



0.3061 1.2743 -0.2630 0.0252 
1.7910 0.3040 
3.6720 
4.0 



directions we use zenith-angle dependence given in Ref. pj[ (see also |20(|"). As the third parame^ 



terization of the atmospheric muon flux we use the formula given in Ref. 
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C. Prompt muon contribution 

Atmospheric prompt muon flux predictions are reviewed in Refs. [3, Q]- Ratios of the differential 
energy spectra of muons at sea level originated from charmed particle decays to that of (vr, K)- 
decays (conventional muons) calculated for a variety of charm production models are shown in 
Fig. [2] (see also (g^). Here PRS stands for the model [sl, GGV for RQPM and QGSM 



for 



18l |. and VZ for Volkova and Zatsepin js^). Among them we dwell below on quark-gluon 



string model (QGSM), as a sample of phenomenological nonperturbative approach, and also on 
some of models based on perturbative QCD computations, GGV [10] and GBW jl4l |. 

Gelmini, Gondolo and Varieschi (GGV) [10] have included NLO corrections for the charm 
production with xg{x) ~ x^'^, (A varying in the range — 0.5). These results obey the following 
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FIG. 2: Ratio of the prompt muon flux to the conventional one at ground level. 



parameterization for the sea- level muon fluxes (see also 22l]): 



(^M^ 'Ml GeV 



- (a+by+cy^ +dy^ ) 



, cm-^s-^sr-^GeV-^ 



(4) 



where y = log^Q(£'^/l GeV). The parameters are given in the Table Ull We choose two represen- 
tative sets corresponding to A = 0.1 (GGVOl) and A = 0.5 (GGV05). 



TABLE II: GGV parameters for the prompt muon fluxes. 



Model 


A, 


— 2 —1 — 1 

cm s sr 


GeV-i 


a 


b 


c 


d 


GGVOl 




3.12 X 10- 


-6 


2.70 


-0.095 


1.49 X 10-2 


-0.2148 X 10-3 


GGV05 




0.58 X IQ- 


-6 


1.84 


0.257 


-4.05 X 10-2 


2.455 X 10-3 



QGSM flux parameterization (that is valid for < 80°) may be written 10\ as 
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E. 



100 TeV 



-3.02 



1 + 



E„ 



100 TeV 



-2.02' 



-0.165 



cm-2s-isr-^GeV~\ 



(5) 



As last representative model, we consider flux calculation within the saturation model proposed 
by Golec-Biernat and WuthsofF [14i |. For this model, we consider two cases (laj): GBWl, where 
the protons are taken to be the primary, and GBW2, where we include the effect of heavy elements 
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also. The sea level prompt muon flux due to GBWl and GBW2 can be parameterized as Eq. 
and Eq. ([7]) respectively: 



rp \ -2.17145-0.0498% 

(E^) = 2.35 • 10-8 ( j , cm-h-hT-'GeY-\ (6) 



/ P \ -1.79371-0.10873J/ 

/,GBW2(^^) ^ 09 . 10-8 f j , cm-2s-isr-iGeV-i 



(7) 



These two cases are different in nature, with the expectation that GBW2 should lead to a decreased 
muon flux at higher energies. 

D. Method to calculate the muon flux under thick layer of the rock 

In these computations we base on the semianalytical method for the solution of muon transport 
equation stated in Ref. [2^ (see also [l^ . [20I). The method allows to consider real atmospheric 
muon spectrum and the energy behavior of discrete energy loss spectra due to radiative and pho- 
tonuclear interactions of muons in matter. Only ionization energy loss of muons are treated as 
continuous one. The method provides effective tool to compute the energy spectra of cosmic-ray 
muons at large depths of homogeneous media. The benefits of this approach are to carry out 
verifications of the primary CR spectrum and composition, charm production models, models of 
the photonuclear interaction with high performance and good precision. This enables to estimate 
the sea- level muon spectrum using the data of underground/underwater measurements evading the 
difficult inverse scattering problem. 

III. EXPECTED MUON FLUX AT THE DEPTH OF PUSHEP SITE 

Zenith-angle distributions of the conventional muon flux calculated for five values of the minimal 
energy of muons in the range 10-10^ GeV at depth 1.3 km of INO detector are shown in Fig. [3l Here 
solid lines represent computations for the surface muon spectrum ^IS] by Bugaev et al. with usage 



of the angle dependence obtained in Ref. [3l (see also 
on solid ones but near horizontal directions, show results 
whereas dotted ones show that for the spectrum by Reyna 



20l|). Dashed lines, almost superimposed 



br the spectrum by Tang et al. 
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2ll |. The geometry of the INO site is 



reflected in the flat shape of the underock distribution (see Fig. [T]). Zenith-angle dependence of 
the conventional and prompt muon fluxes at the INO depth are shown in Figs. H] and [5l In Fig. H] 
are shown the prompt muon flux at muon energy above 100 and 200 TeV calculated with QGSM 
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FIG. 3: Angle distributions of the conventional muons near the INO detector. 
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FIG. 4: Zenith-angle distributions of atmospheric muons at 100 and 200 TeV. 

charm production cross sections (dash-dotted lines) and that of GGV models. For muon energy 
above 500 TeV we also plot predictions obtained for GBW model (dashes line in Fig. [5|). As one 
may clearly observe in Fig. [5l measurement of high muon flux near the vertical at INO depth could 
allow to discriminate between GGV01(A = 0.1) model and GGV05(A = 0.5) or QGSM one. While 
the GBW prompt muon flux is unlikely to be observed at 500 TeV. 

Differential muon spectra (left panel) at INO depth and integral ones (right panel) are presented 
in Fig. [6] for cos^ = 0.7, where solid line shows the conventional muon flux obtained with usage of 
Bugaev et al. boundary spectrum and circles denote that for Gaisser's spectrum. We can see in 
Fig. [6] that crossover energy for the conventional muon flux and the GGV05 prompt one is about 
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300 TeV, therefore it seems that more suitable for the prompt muon identification is to analyse 
the zenith-angle dependence of high-energy muon flux (see Fig. S]) . 
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FIG. 5: Very high-energy zenith-angle distributions of atmospheric muons. 
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FIG. 6: Energy spectra of the atmospheric muons near the INO detector. 



Number of the muon events per steradian per year expected at INO detector near direction 
cos^ = 0.7 is presented in Table Hill (see details in Ref. |12l|). Last three columns in Table Hill 
represent the ratio of the conventional muon flux to the prompt muon one due to three charm 
production models, GGVOl, GGV05 and QGSM, respectively. 



IV. SUMMARY 



The shape of zenith-angle distributions of conventional muons is nearly flat (see Figs. [SHS]). 
Therefore muons arriving at the detector close to vertical directions are more favorable to measure 
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TABLE III: Number of the muon events per steradian per year expected at INO detector. 



E^, TeV 


conv. 


GGVOl 


GGV05 


QGSM 


^GGVOl ^GGV05 


^QGSM 


10 


60097 


1235 


1353 


3037 


0.02 


0.022 


0.05 


50 


832 


73 


105 


159 


0.087 


0.126 


0.19 


100 


132 


20 


34 


41 


0.15 


0.258 


0.31 


200 


20 


5.0 


10 


10 


0.25 


0.50 


0.50 


300 


6.0 


2.0 


5.0 


4.0 


0.33 


0.83 


0.66 


400 


2.6 


1.0 


2.6 


2.2 


0.38 


1.0 


0.85 


500 


1.4 


0.6 


1.6 


1.2 


0.43 


1.14 


0.86 



the prompt muon flux. The prompt muon contribution to the atmospheric muon flux increases 
with energy because of lower value of the energy spectrum index. The "crossover" energy, Ec, 
at which the prompt muon Hux becomes equal to the conventional one, depends strongly on the 
charm production model. Following numbers can illustrate (see Fig. [2]) the at INO depth for 
some of charm hadroproduction models: .B^^^os ^ 250 TeV, ^Q^^^ ~ 300 TeV, ^^^^oi ^ qqq 
TeV. 

From the Table Hill we can see that prompt muon flux contribution due to GGVOl model, for 
example, may differs from that for the GGV051 model (or QGSM) by factor 2 at E"^ > 200 TeV. 
In other words, expected number of muon events inside the INO detector may increase by 50 % at 
the energy above 200 TeV if GGV05 or QGSM predictions are reasonable. 
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